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Abstract 

This work aims at illustrating that, for a class of leading- twist hard exclusive reactions involving 
two heavy quarkonia, the light-cone approach, when equipped with the strategy of refactorization 
of the light-cone distribution amplitude of quarkonium, can be employed to elegantly reproduce 
the corresponding predictions made in the nonrelativistic QCD (NRQCD) factorization approach, 
order by order in perturbative expansion. Taking the electromagnetic form factor of the Be meson 
at large momentum transfer, , as a concrete example, we compare the results obtained from both 
NRQCD-based and light-cone-based calculations through the next-to- leading order (NLO) in q^, 
while at the leading order (LO) in both velocity and expansion, and explicitly confirm their 

mutual agreement. As a byproduct, we apply our NLO result to explore certain features about the 
asymptotic behavior of the heavy-light meson form factor. We also address the major theoretical 
obstacles that prevent us from establishing an analogous equivalence between these two approaches 
for the double charmonium production process of phenomenological interest, e+e" — >• J /if) + rjc- 
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I. INTRODUCTION 



The utility of perturbative QCD crucially rests upon the idea of factorization. For a typi- 
cal QCD process with large momentum transfer (either inclusive or exclusive) , factorization 
provides an essential tool to systematically separate the short-distance, perturbatively cal- 
culable effects from the long-distance, yet universal, nonperturbative effects. For a hard 
exclusive processes involving a few hadrons, the well-known coUinear factorization (also re- 
ferred to as light-cone approach in literature) [1, 2], fulfills such a separation by expressing 
the amplitude as the convolution of the hard-scattering part with the nonperturbative yet 
universal light-cone distribution amphtudes (LCDAs) of the corresponding hadrons. The 
classic applications of light-cone factorization are exemplified by the 71 — 7 transition form 
factor and tt electromagnetic (EM) form factor [1, 2], and the nonleptonic B meson de- 
cays [3, 4], to which a vast amount of literature has been devoted. 

Recent advancement of the high-luminosity high-energy collider facilities renders the 
study of hard exclusive processes involving heavy quarkonium a fertile research frontier. 
Perhaps a great amount of interests toward this topic have been triggered by the obser- 
vation of several double-charmonium production processes in two B factories several years 
ago [5, 6]. 

For a hard exclusive process comprising entirely of light hadrons, such as pion EM form 
factor, only the hard momentum transfer Q and the nonperturbative QCD scale, Aqcd, 
play dynamical roles. By contrast, hard exclusive reactions involving quarkonium bring 
forth much more complexity as a number of new dynamic scales have come into play alto- 
gether. Among the important scales inherent to a heavy quarkonium, are the heavy quark 
mass, ni, the typical momentum of heavy quark, mv, and the typical binding energy of a 
quarkonium, mv^, where v denotes the typical velocity of heavy quark inside a quarkonium. 
As a consequence, it is theoretically more challenging to analyze the hard reactions involving 
heavy quarkonia than those involving only light hadrons. 

Provided that the characteristic momentum transfer, Q, is much greater than m, the 
light-cone approach undoubtedly remains to be a valid and powerful method, since it is a 
formalism based on the expansion in powers of 1/Q. Nevertheless, there also exists an- 
other influential factorization approach, the NRQCD factorization [7], which combines the 
effective-field-theory (EFT) machinery and the hard-scattering factorization. Tailor-made 
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to tackle quarkonium production and decay processes, this method exphcitly exploits the 
nonrelativistic nature of quarkonium by factoring the perturbative quantum fluctuations of 
distance 1/m or shorter, from the nonperturbative effects governing the transition of a heavy 
quark pair into a physical quarkonium, which occurs at a distance of 1/mv or longer. As 
a consequence, NRQCD factorization allows one to express the amplitude of an exclusive 
quarkonium production reaction in terms of an infinite sum of products of short-distance 
coefficients and the vacuum-to-quarkonium NRQCD matrix elements, whose importance has 
been organized by the powers of v. 

Both factorization frameworks, despite spothghting different dynamic aspects, are well- 
established approaches derived from the first principles of QCD. They arc commonly viewed 
as two drastically different, and, even repulsive methods. For example, in recent years, 
intensive investigations on the reaction e'^e~ J/ip + rjc have been conducted from both 
frameworks [8-15], and there have been some disputes on which approach is superior when 
applied to the hard exclusive reaction involving quarkonium [14, 16]. 

However, it seems fair to state that, when coping with hard exclusive process involving 
quarkonium, each approach has its own strength and weakness. For example, in NRQCD 
factorization, the short-distance coefficients at each order in would inevitably contain large 
logarithms of the ratio of two vastly different scales, Q {e.g., the center-of-mass energy for 
quarkonium-|-7 or double-charmonium production in e"*"e~ annihilation) and m. To improve 
the stability of perturbative expansion, a systematic disentanglement of these two scales is 
clearly called for. On the other hand, the light-cone approach, being a general formalism 
applicable to any type of hadrons, usually does not sufficiently exploit the nonrelativistic 
trait of quarkonium. That is, the LCDAs of quarkonium are often determined on purely 
phenomenological ground [14, 16-18]. The unsatisfactory feature of this approach is that, 
quarkonium LCDAs inevitably encompass the coUinear degrees of freedom of vastly varying 
virtualities, some of which may be high enough to be perturbatively calculable and separable 
from the remaining parts. 

Recently it has been reahzed that these two methods need not be regarded solely as 
rivals, rather they can be tied coherently to compensate the shortcoming of each other. 
In Ref. [19], it is suggested that the light-cone approach can be utilized to refactorize the 
NRQCD short-distance coefficients for a class of single-quarkonium production processes, 
exemplified by 775— 7 form factor and Higgs bosn decay to T -|- 7, by which the scales Q and 
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m can get fully disentangled. Consequently, the coUinear logarithms of type Q;"ln"((5^/m^) 
can be readily identified and resummed in this framework ^. Furthermore, it is envisaged in 
[19] that the hght-cone approach can serve as an alternative means to reproduce the NRQCD 
factorization predictions at the leading order in 1/Q and order by order in a^. Therefore 
this refactorization procedure can spare a great amount of labors compared with the direct 
higher-order computations in NRQCD framework. 

There has already existed another interesting development from a different but com- 
plimentary perspective, which explores the possibility of further factoring the LCDAs of 
quarkonium [21]. The underlying idea is that, since the virtualities of the coUinear modes 
encoded in quarkonium LCDA can range from 0{m^) down to 0{m^v'^) or lower, it seems 
profitable to identify and separate the shorter-distance coUinear quantum fluctuations out 
of the quarkonium LCDA. Thus, it is suggested in [21] that, the LCDAs of a quarkonium 
can be matched onto an infinite sum of the product of the perturbatively-calculable univer- 
sal jet functions and the vacuum-to-quarkonium NRQCD matrix elements, organized by v 
expansion. At the lowest-order (LO) in the corresponding jet functions associated with 
the twist-2 LCDAs of the S-wave quarkonia, such as r^c, and Sc, have been calculated 
through the next-to-leading-order (NLO) in [21, 22]. 

Both versions of refactorization, operating either on NRQCD short-distance coeffi- 
cients [19], or on quarkonium LCDAs [21, 22], aim to achieve an optimized scale separation 
through bridging the hght-cone and NRQCD approaches together. Though motivated from 
different considerations in curing the particular shortcoming of the respective factorization 
approach, both of them in fact lead to equivalent final predictions, that is, the exclusive re- 
action amplitude will be expressible in the form of a hard scattering kernel convoluted with 
the pertubative jet function, then multiplied by the nonperturbative vacuum-to-quarkonium 
NRQCD matrix element. Practically speaking, the second version of refactorization strat- 
egy [21], which highhghts the universality of the perturbatively-calculable jet functions, is 
presumably more convenient to employ, and we will follow this perspective in this work. 

The main motif of this work is to verify the correctness and effectiveness of the refactor- 
ization strategy for a class of leading-twist hard exclusive processes involving two quarkonia, 

^ It is interesting to remark that, resummation of the large kinematic logarithms for the above-mentioned 
exclusive single-quarkonium production processes have also been performed by Shifman and Vysotsky 
thirty years ago [20], long before the inception of the NRQCD factorization approach. 
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which is more nontrivial and interesting than those single-quarkonium production processes 
considered in [19] . Specifically, by employing the NRQCD approach and the hght-cone ap- 
proach separately, we will investigate the EM form factor of the meson at asymptotically 
large momentum transfer. In spite of lacking urgent phenomenological incentive, this process 
can serve as an ideal theoretical laboratory to corroborate our understanding in a nontrivial 
way We will explicitly verify that, to the accuracy at LO in v but through NLO in as, the 
light-cone approach, when equipped with the machinery of refactorization, can be utilized 
in an effective and elegant manner to reproduce the asymptotic NRQCD predictions for the 
Be EM form factor. Aside from the calculational advantage, the refactorization strategy 
may also shed light on how to ameliorate the severe scale dependence associated with the 
NLO NRQCD predictions observed in various quarkonium production processes. 

Apart from its great utility, however, it is worth emphasizing one important caveat for 
this refactorization program. Due to some long-standing problems in the coUinear fac- 
torization framework, the refactorization strategy, at its present formulation, can not be 
applied successfully to any types of hard exclusive reactions involving quarkonium. As we 
will see later, the phenomenologically interesting double-charmonium production process, 
e^e^ J/ip + r]c, constitutes such a very example. 

The rest of the paper is structured as follows. In Sec. II we define the electromagnetic 
form factor of the Be meson. In Sec. Ill, we compute the Be EM form factor through the 
NLO in as while at the LO in velocity within NRQCD factorization framework. In Sec. IV, 
we reinvestigate the Be form factor at large in the framework of coUinear factorization 
that implements the refactorization strategy, to the leading power in l/Q'^ and v^, yet 
through the NLO in a^. We then compare the light- cone- based prediction and the asymptotic 
NRQCD result through the NLO in as, and explicitly establish their equivalence at large 
Q^. In Sec. V, we apply our knowledge of the NLO expression for the Be form factor 
to discuss certain features of the EM form factor of a heavy-hght meson like and 
speculate on possible refactorization procedure suitable to this situation. In Sec. VI, we 
discuss the pecuhar asymptotic behavior of the NLO NRQCD prediction to the reaction 
e~^e~ J/ip + rje, and briefly address the inapplicability of the refactorization strategy to 
this case. Finally we summarize in Sec. VII. 
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FIG. 1: Schematic depiction of the quark transition process relevant to the Be EM form factor. 

II. DEFINITION OF THE ELECTROMAGNETIC FORM FACTOR 

Imagine a Be meson struck by an electromagnetic probe undergoes an elastic scattering. 
The information of this elastic scattering is encoded in the EM form factor of B^, dubbed 
FiQ"^) in this work. It can be introduced in full analogy with the EM form factor of the 
charged pion: 

{Bnp')\Ji:jBt{P)) = F{Q^){p+pr. (1) 

where the electromagnetic current = Yli^q'il'^'i (^g = +| up- type quarks and = 
— I for down-type quarks), and the momentum carried by the EM probe is q = P' — P, 
with = —q^ > signifying the momentum transfer. The structure of (1) is uniquely 
determined by Lorentz invariance and the EM current conservation. A cartoon that depicts 
the Be form factor is sketched in Fig. 1. 

Throughout this work, we will be primarily interested in the asymptotic behavior of the 
form factor in the limit of ^ M^^ ^. In such a situation, both NRQCD factorization and 
coUinear factorization are valid theoretical tools to predict this form factor. Although the 
dimensionless form factor itself is a Lorentz scalar, the physical picture is probably more 
transparent if the Breit frame is specifically borne in mind. In such a frame, a Be meson 
moving very fast along the z axis gets hit by a highly virtual photon, then heads back along 
the —z axis with equal speed. 

^ At present there are no experimental facilities which can directly measure this form factor at large spacelike 
momentum transfer. The next generation e~^e~ machines, e.g. the prospective International Linear 
Collider, or the super-high-luminosity machine operating at the Z'~' pole, may be possible to access this 
form factor in the time-like region. 



III. EM FORM FACTOR OF IN NRQCD FACTORIZATION 



In a hard exclusive process involving quarkonium, a pair of heavy quark and an antiquark 
have to be created/annihilated in short distance; and in order to have a significant probabil- 
ity to form/disintegrate a quarkonium, the typical relative motion between the quark and 
antiquark should be necessarily slow. The first condition guarantees that the asymptotic 
freedom can be invoked to compute the hard-scattering quark amplitude in perturbation 
theory. The second condition implies that, this hard-scattering amplitude is insensitive to 
small change of the relative momentum of the pair, thus one may expand the amplitude in 
power series of v, subsequently absorb those I'-dependent factors and the quarkonium wave 
function into the nonperturbative matrix elements, which are dictated by the quarkonium 
binding mechanism. NRQCD factorization then naturally arises from this physical picture. 

At the LO in v, the EM form factor in NRQCD factorization can be written in the 
following form: 



i^NRQCD(g') = C{Q;m,,m,) | -^j + 0{v') , (2) 

where C{Q\mc^mij) denotes the corresponding short-distance coefficient. /^^J represents the 
nonperturbative vacuum-to- 5c NRQCD matrix element, which is defined through 



f(0) 
JL 



AO) ^ 

JBc — 



{0\xli^c\B^) ^ m\4xb\0)\ = ]I^^bM- (3) 



Me 

Here A^c = 3 is the number of the colors, ipc and Xb denote the Pauli spinor fields associated 
with the c and b fiavors, respectively. The Be state in the above NRQCD matrix elements 
is normalized non-relativistically. As an alert reader may readily tell, f^^ in fact is nothing 
but the Be decay constant rephrased in the NRQCD context ^. The superscript reminds 
that this identification is accurate only at the lowest order in ag and v. As displayed in (3), 
this entity can also be linked with ipBci^), familiar Schrodinger wave function at the 
origin for Be- Furthermore, we will assume Mb^ — mc + rrib everywhere in this work, which 
is legitimate to the LO accuracy in v expansion. 

As stressed before, we will concentrate on the kinematic situation Q ^ rric^i, ^ Aqcd- 
The NRQCD short-distance coefficient, C{Q; rric, rrib), can then be deduced reliably through 



^ Note that our convention of defining the decay constant of a charged pseudoscalar meson is such that 
U = 132 MeV. 



7 



the perturbative matching procedure. Consequently, it can be expanded in power series of 

as C = C(°) + ^C(^) H . Accordingly, the NRQCD prediction to the form factor 

can be organized as 

p _ p(0) O'silJ'R) p(l) . - X 

-f^NRQCD - -f'NRQCD H "^NRQCD " l^j 



A. NRQCD prediction at LO in ag 

At tree level, there arc only four Feynman diagrams for the quark-transition process 
7* + cb{^Si^\P) cbCsi,^\P'). They can be obtained by replacing the shaded ellipse in 
Fig. 1 with a single gluonic exchange between c and b quarks, as well as by intersecting the 
EM current either to the c or 6 quark line. Assuming that both constituents in each pair 
have vanishing relative motion, it then is a straightforward exercise to deduce (7^°^: 

>{Q; me, rrib) = 47rCj.Q;,(//^) le^ ^ - ^^^^^3 1 . (5) 

For abbreviation, we have defined Xq = and xo = 1 — xq — ^. Cp — denotes 
the Casmir of the fundamental representation for SU{N^ group. 

In order to make contact with the leading-twist light-cone predictions, which will be 
reported in next section, we are particularly interested in the asymptotic behavior of Eq. (5) 
in the limit Q > Mb, ^: 



C•<°>„(«;m„m.) = ^^^f^^(|-|). 



(6) 



The NRQCD prediction correctly embraces the l/Q^ scahng of the form factor of a pseu- 
doscalar meson, as it should be. 

Without sacrifice of the essential goal, and for great technical simplicity, it is also useful to 
imagine a gedanken experiment by playing with a fictitious Be meson, whose both constitutes 
carry equal mass: rric — nib. In such an idealized case, xq — xq — ^. Eq. (5) then reduces 



Throughout this paper wc adopt the notation that x = 1 — x for any momentum fraction x. 
^ In the case for the physical meson, we will not pursue the further scale separation between rric and 
nib, assuming no significant hierarchy between xq and xo- Nevertheless in Section V, we will explore the 
consequence of the heavy-quark limit mt '> nic, which might be relevant for the B meson form factor. 
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to 6 

C(°)(g;me,me) = 167rC^e^+a,(/x^)^^-^^, (7) 

where 6^+ = Cc — = +1 is the electric charge of the B'^ meson. Its hmiting behavior, 
C'i^ym((5; mc) = I&'kCfc^s/Q^-i is particularly simple. 

B. NRQCD prediction at NLO in as 

To assess the NLO perturbative correction to the NRQCD short-distance coefficient, C^^\ 
one can follow the standard matching procedure by computing the on-shell quark amplitude 
for 7* + c6 — )> cb to NLO in a^, neglecting the relative momentum between c and h for both 
incoming and outgoing Be- Excluding counter-term diagrams, there are in total 64 NLO 
Feynman diagrams. Since several distinct scales: Q, mi,, rric, will enter in loop integrals, it 
is conceivable that fulfilling all the analytic management would be a highly formidable task. 

In this work, we will utilize one of the world-leading automated Feynman Diagram 
Calculation packages (FDC), to fulfill all the required algebras of tensor-reduction and 
one-loop scalar integral. FDC is a powerful program based on the symbolic language RE- 
DUCE, which was designed to automate the perturbative quantum field theory calculation 
in computer. FDC was initially developed by one of us long ago [23], and the function of 
automatic one-loop calculation has recently been reahzed by Gong and Wang [24] . In recent 
years, the FDC package has been vividly applied to numerous quarkonium production and 
decay processes, and has withstood many highly nontrivial tests [12, 25-28]. 

The masses of the c and b quarks have been retained explicitly in the calculation. The 
nif = 3 species of light quark flavors: u, d, s, which occurs in the gluon vacuum polarization 
diagrams, have been treated as massless. Because the masses of the b and c are still much 
smaller than the momentum transfer Q, one should count the number of active flavors as 
Uf ^ riif + 2 ^ 5. 

Both ultraviolet (UV) and infrared (IR) (in our case including soft and Coulomb) diver- 
gences may appear in an individual diagram, while the would-be coUinear divergence has 

^ Even with rub set equal to rric, such a fictitious Be state, as a flavor non-singlet, should be distinguished 
from the flavor-singlet pseudoscalar quarkonium such as rjc- The EM transition 'y*r]c rjc is simply 
forbidden due to violation of the C-parity. 
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been cutoff by tfie lieavy quark mass. Dimensional regularization fias been used to regularize 
both UV and IR divergences ^, and the relative velocity v used to regularize the Coulomb 
singularity. All the occurring UV divergences are canceled systematically by incorporating 
the counterterm diagrams. As usual we adopt the MS scheme to renormalize the QCD cou- 
pling constant. One has the freedom to choose some specific renormalization prescriptions 
for the quark and gluon fields, e.g., MS or on-shell scheme. However, the LSZ reduction 
formula guarantees that the resulting short-distance coefficient C^^\ which is inferred from 
the on-shell quark amplitude, is free from any scheme ambiguity, at least to the NLO which 
we are considering. 

It turns out that the soft IR divergences cancel upon summing all the diagrams, and, as 
usual, the remaining Coulomb divergences can be factored into the NRQCD matrix element 
through the matching. Hence the final expression of C-^-* becomes both UV and IR finite, 
depending only on Q, rric, rrib, and the renormalization scale, /ir, respectively. 

In the case of a physical Be state with rric ^ mi,, this calculation poses a rather nontrivial 
challenge to the capability of the FDC package, which is perhaps by far the most involved 
NLO calculation for the hard exclusive QCD process. After improving the current algorithm 
for the scalar-integral part, and passing several nontrivial internal checks, FDC can indeed 
successfully generate the correct result. Unfortunately, the resulting analytic expression for 
C^^^ turns to be extremely involved and pathologically lengthy. There seems absolutely 
no way to directly manipulate on this analytic result, not even mentioning to deduce its 
asymptotic behavior. Therefore we must be contented with providing only the numerical 
results. 

In the left panel of Fig. 2, we explicitly show the ratio of the NLO form factor to the LO 
form factor of a physical Be, -^nrqcdZ-^nrqcd (which equals to C'^'^^/C^^^) in a wide range 
of Q. For simplicity we have set /ir — Q, to eliminate the potentially large UV logarithm. 

The most noteworthy feature is that, after Q > 400 GeV, the ratio -^nrqcdZ-^nrqcd starts 
to grow linearly with InQ (or more precisely, with \n{Q/MBj to balance the dimension). 
This particular type of single-logarithm scaling is the very vestige of those would-be coUinear 
divergences that had been cutoff by the nonzero b and c masses. Within the confine of 



During the calculation we need to handle the trace of a string of Dirac Gamma matrices containing two 75 
in Z) = 4 — 2e dimension. As pointed out in [29], one can safely utilize the naive dimensional regularization 
prescription to annihilate these two 75. 
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10 10 10 10 10 10 10 10 10 10 

Q (GeV) Q (GeV) 

FIG. 2: The ratio F^j^^ (Q^) / F^^ {Q^) as a function of Q with Mb, = 6.3 GeV, n/=5 (/3o = f ), and 



Aifl = Q- Both the NRQCD and hght-cone predictions are shown, where the former is represented 
by sohd line, and the latter by the dashed line. The left panel is for the EM form factor of the 
physical Be state with nic = 1.5 GeV and = 4.8 GeV, while the right panel for a fictitious Be 
state with rric = ruf, = 3.15 GeV. Numerically, the dashed line in the left panel can be parameterized 
by 0.680 In Q + 1.624, and that in the right panel by 2.152 In Q — 1.795, where Q is in the unit of 
GeV. 

NRQCD approach, it is difficult to ascertain the explicit form of this single logarithm without 
a complete NLO calculation. Nevertheless, it has been illustrated that [19], the light-cone 
approach, can be used effectively to identify, and resum, such leading collinear logarithms 
to all orders in for a class of single- quarkonium production processes. In the following 
section, we will explicitly see how the refactorization strategy can help to reproduce this 
logarithm in a rather straightforward manner. 

In order to acquire a clear view on the asymptotic behavior of C^^\ it is helpful to examine 
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the NLO QCD corrections to the EM form factor of the aforementioned fictitious Be state. 
The analytic expression of (7*^^^ in this case, though still involved to some extent, fortunately 
becomes far simpler and much more manageable than that in the case of the physical B^.. We 
note that the complexity of this calculation is comparable with that of the NLO correction to 
J/ip + rjc EM form factor. To our purpose, it is desirable to deduce the asymptotic behavior 
of C^^^ in a closed form. After some judicious yet tedious manipulations on the analytic 
output of FDC, we find the asymptotic expression of C^^^ for this fictitious Be state to be 

Here /3o — ^-Nc — ^Uf is the one- loop coefficient of the QCD /3 function. The occurrence 
of the /3o In hr term in (8) is as expected, required by the /Xfl-independence of the NRQCD 
short-distance coefficients. 

In the right panel of Fig. 2, we have shown that the ratio -^nrqcd/ -^nrqcd fictitious 
Be state over a wide range of Q. The exact NLO prediction is juxtaposed together with the 
asymptotic expression, (8). As indicated clearly in Fig. 2, as long as Q > 100 GeV, this 
asymptotic expression starts to converge to the exact NLO result quite well. 

Note the coefficient of the coUinear logarithm in ci^ym in (8) is of the algebraic structure 
oc Cf(3 — 2 In 2), similar to what appears in the 77c — 7 transition form factor [19, 20, 30]. 
This will be easily understood in the light-cone-based framework in next section. 



IV. EM FORM FACTOR OF Be IN LIGHT-CONE APPROACH 

At asymptotically large Q^, both the incoming and outgoing Be, when viewed in the Breit 
frame, move nearly with the speed of light. By virtue of the asymptotic freedom of QCD, 
the hard-scattering quark amplitude can be accessed by perturbation theory. Since both the 
quark and antiquark inside Be are dictated by the light-like kinematics, the hard-scattering 
amplitude is insensitive to the small change of the quark mass, nicb, as well as the transverse 
momentum carried by the quark and antiquark, p±. Thus, the amplitude can be expanded 
in powers of p± and me,b, while the nonperturbative wave function together with the p± and 
quark-mass-dependent effects can be lumped into the LCDAs of Be- This picture naturally 
endorses the applicability of coUinear factorization. 

At the leading power in l/Q, the EM form factor of Be can be factored into the convolution 
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of the perturbatively calculable hard-scattering amplitude Th with the leading-twist LCDAs 
of Be, signified by ^bS^)- 

FLc(g') = ['dx['dy^*j,^{y,f,l)TH{x,y,Q',fil,fil)^BAx,f,l) + 0{l/Q'), (9) 
Jo Jo 

where x, y represent the fractions of light-cone momentum carried by the c quark in the 
incident and outgoing Be states, and /ir, /ip denote the renormalization and factorization 
scales, respectively. 



A. Outline of refactorization strategy in light-cone framework 

The factorization theorem (9) warrants that, the hard-scattering kernel, Th, can be sys- 
tematically improved in perturbation theory. When computing the hard-scattering ampli- 
tude for Y + c{xP)b{xP) c{yP')h{yP') to the leading-power accuracy, the masses and the 
transverse momenta of the quarks and antiquarks have been neglected, so both c and h move 
parallel to the momentum of Bf.. Obviously, the Th in our case would be the exactly same 
as the corresponding one in the tt form factor, which has been available for a long while. It 
is convenient to organize the hard-scattering amplitude in power series of a^: 

TH{x,y,Q\i,l,i,l) = TP{x,y,Q') + ^^^T^^\x,y,Q\^il,^l)-r---. (10) 

TT 

A crucial ingredient of the factorization formula (9) is the nonperturbative LCDA of the 
Be meson. The Be LCDA are conventionally parameterized as follows: 

^BA^,^J^l) = ^^J{x,f,l), (11) 

where /b^ is the (physical) decay constant of Be, and (f) can be viewed as a function character- 
izing the probability amplitude for the c (b) quark to carry the fractional light-cone momen- 
tum X (x). The factorization scale /ip also enters into 0, to preclude those shorter-distance 
configurations with transverse momentum greater than fip to be erroneously included in 
the nonperturbative LCDA. The local limit of the LCDA imposes a model-independent 
normalization condition j^dx(j){x, ^'^p) = 1, valid for any ^p. 

In the case of pion, the distribution amplitude 0, which necessarily probes the long- 
distance fluctuation of order Aqcd, is a genuinely nonperturbative object. Therefore it can 
only be tackled by nonperturbative tools such as lattice QCD simulation [31, 32], or by 
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phenomenological methods such as the QCD sum rules [33, 34], or the holographic QCD 
models [35] that have recently gained much popularity. Unfortunately, despite many inten- 
sive studies for decades, the pion distribution amplitude is still not accurately known even 
today. 

By contrast, according to the philosophy of refactorization, the distribution amplitude 
for a quarkonium may be viewed as the short-distance coefficient (jet) function associated 
with matching the quarkonium LCDA onto the NRQCD vaccum-to-quarkonium matrix 
element, while the nonperturbative effect of order mv or lower has been entirely encoded 
in /b^. Since heavy quark mass can serve as an infrared cutoff, this jet function, which 
encompasses the effects of collinear modes with virtuality of order m^j, ^ ^qcd> should be 
reliably accessible by perturbation theory owing to asymptotic freedom. At the LO in v, 
this function can be expanded in powers of as'. 

4>{x,n%) = 0(0) (a:) + (a:, /.!) + •••. (12) 

TT 

The LO jet function can be trivially inferred, 

0(°)(x)=(5(x-xo), (13) 

which simply reflects that both c and b partition the momentum of Be commensurate to 
their mass ratios. This is compatible with the LO NRQCD expansion, in that c and b are 
at rest relative to each other in the Be rest frame ^. 

Separating the nonperturbative LCDA of a meson as the product of the decay constant 
and a distribution function i> in (11), is merely a parametrization in the case of pion. 
Nevertheless, such a parametrization for a quarkonium LCDA, which is referred to as refac- 
torization in this work, is a rather useful concept and can lead to quite nontrivial outcome. 

In line with the spirit of refactorization, all the nonperturbative aspects of the twist- 
2 LCDA of Be, is evidently contained in the Be decay constant /s^. In order to make 
contact with the NRQCD prediction in Sec. Ill, it is useful to specify the connection between 
the physical Be decay constant and f^^ introduced in (3). Beyond the LO in a^, it is 
necessary to distinguish fs^ from its approximate zeroth-order result, Their connection 

^ Note that the LCDA of Be is asymmetric about .t = i, which is rather different from the n case. The 
reason can be easily traced, that is because the flavor symmetry between heavy quarks c and 6, unlike the 
isospin symmetry between u and d, is badly broken. 
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is best framed in the NRQCD context, through matching the QCD axial- vector current to 
its NRQCD counterpart order by order in ag and in v. Schematically, the conversion relation 
between fs^ and can be expressed as follows: 

/b. =/g'(l+^^^fS + ■■■), (14) 

which is valid only at LO in v. For such a static quantity, it is customary and most appro- 
priate to choose the renormalization scale for ag appearing in (14) to be around M^^, which 
characterizes the typical virtuality of quantum fluctuations integrated out by NRQCD. 

Analogous to (4) in NRQCD factorization, we also organize the light-cone prediction to 
the form factor of B^. in the power series of ct^, 

F.c = FS + ^fS + ---. (15) 

TT 

Once the Th, and /^^ are separately known through the NLO accuracy in a^, following 
(9), we can then readily recognize the LO and NLO predictions to the Be form factor in 
light-cone perturbation theory. Obviously, the LO form factor simply assumes the following 
form: 

F^g ^ 0(°)®r|')®0(°), (16) 

where (8> designates the convolution, and for simplicity we have suppressed the common 
multiplicative factors, such as the square of f^^\ 

The NLO contribution to the form factor, i^L^, would receive contributions from several 
different sources: 

F^^^ ~ 0(°)®rg)®0(°), 0W®4°)®<^w, fg0(°)®rf 

(17) 

In (15), we have deliberately kept some ambiguity on the choice of the strong coupling 
constant, that is, we have not specified the scale at which the ag should be evaluated. It is 
worth noting, however, the different that stem from different NLO sources in (15), are 
in principle affiliated with quite different scales, e.(/., fip, and Mb^, as clearly indicated 
in (10), (12) and (14). Hence there seems no rationale to attach a uniform scale to all the 
occurring ccg, contrary to the form suggested by (15). 
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Since the primary goal in this work is to make a critical comparison between the light-cone 
and NRQCD approach at the NLO level, it seems tolerable, at this stage, not to meticulously 
distinguish the scale attached to each a^. For the sake of making contact with the NRQCD 
prediction, and also for simplicity, from now on we will assume all the in (15) to be 
evaluated at a single scale, say, /ir. The error induced by such a mismatch will propagate 
to the higher order. 

At the end of this section (Sec. IV E), we will envisage a more appropriate scale setting 
scheme, and briefly address the possible advantage of the refactorization strategy over the 
conventional NRQCD factorization calculation in ameliorating the scale dependence for a 
fixed-order calculation. 



B. Light-cone prediction at LO in ag 

At tree level, there are four Feynman diagrams, which can be obtained by replacing the 
shaded ellipse in Fig. 1 with a single gluon line between c and b quarks, as well as by attaching 
the EM current either to the c or the b quark line. Assuming the incident and outgoing c 
quark to move coUinear to the respective Be mesons, and taking their corresponding light- 
cone momentum fractions to be x and y, respectively, it is a straightforward exercise to get 
the Born-order hard-scattering kernel: 



\xy xyj 

This expression is known ever since the pioneering work by Brodsky and Lepage [1] . 

As has been pointed out, the LO jet function associated with the Be LCDA in (11) simply 
is (j)^^\x) = 6{x — xq), compatible with the NRQCD expansion at the zeroth-ordcr in v. 

Inserting the explicit expressions of (f)^^\x) and the Born- level hard-scattering kernel, 
(18), into (16), we obtain the LO prediction to the Be form factor: 

Not surprisingly, one finds the exact agreement between this result and the asymptotic 
expression of the LO NRQCD prediction in (6). 
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C. Light-cone prediction at NLO in as 



One may feel that, the agreement between NRQCD and hght-cone predictions at tree level 
is easy to envision, and, more or less trivial. However, concerning the frightening complexity 
of the NLO calculation in NRQCD side, one may agree that, if the equivalence between these 
two approaches at NLO can be further established, it would be an unambiguous indicator 
that the strategy of refactorization is on the correct track. 

According to the schematic recipe (17), one can identify those different components of 
NLO corrections concretely and assemble them together. Consequently, the NLO hght-cone 
prediction to Be form factor, i^LC' ^® expressed as 

- (cc Cfe, xo ^ xo) I . (20) 

T^"* is related to the NLO hard-scattering kernel whose precise meaning will be specified 
in (21). {x~^)^^^ signifies the first inverse moment of the light-cone momentum fraction of 
the b quark, where the superscript imphes that 1/x needs to be folded with the NLO jet 
function. Its concrete expression will be presented in (27). 

To complete the NLO analysis in light-cone framework, we need to work out all the 
involved ingredients in (20) successively. 



Nc \x' 



1. NLO correction to the hard- scattering kernel: T^' 

One indispensable ingredient for a complete NLO analysis is the NLO hard-scattering 
kernel, As was mentioned, this quantity is exactly identical to the corresponding one 
for pion EM form factor, so there is no need to recalculate it in this work. 

The NLO correction to the hard-scattering amplitude for the pion EM form factor has 
been investigated by many authors over the last three decades, but the history of this 
study seems to have followed a somewhat twisted path. The NLO calculation was initially 
carried out by three independent groups in early 1980s. But unfortunately, these results 
did not fully agree with each other [29, 36, 37]. In 1987, after scrutinizing the previous 
calculations, Braaten and Tse were able to locate the origin of the discrepancies among 
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the earlier works [38] . One decade later, after critically reexamining all the existing works, 
Melic, Nizic, and Passek then presented an ultimately consistent version, and also conducted 
a comprehensive phenomenological study by including the evolution effect of pion LCDA [39]. 
In this work, we will take the expression of T^^ from Ref. [39]. 

Prior to quoting the concrete expression, we would like to recapitulate some noteworthy 
aspects encountered in this NLO calculation. Each individual NLO diagram may contain 
the single UV pole, as well as the single or double IR pole. The occurring UV divergences 
can be eliminated by the standard renormalization procedure. Conventionally, both the field 
strength and the strong coupling constant are renormalized within MS scheme. The single 
IR pole may be of soft or collinear origin, and the double IR pole stems from the overlap 
between soft and collinear singularities. Upon summing up all the NLO diagrams, the double 
IR poles and single soft IR poles cancel, while only the single collinear IR pole survives in 
the final expression of the NLO quark amplitude ^. According to the coUinear factorization 
theorem, a standard prescription is to employ the MS factorization scheme to absorb this 
collinear singularity into the pion LCDA, consequently the ultimate expression for the hard- 
scattering kernel, becomes both UV and IR finite. T^^ can be parameterized as 

(21) 

where represents the reduced hard-scattering kernel with EM current attached either 
to the c or to the b. Note T^'' is a dimensionless quantity. The explicit dependence of 
on nn and hf embodies the vestige of those original UV and collinear IR singularities. 

To the intended accuracy, it is sufficient for T^^ to convolve with two LO LCDAs 4>^^^ of 
Be, which are (5-functions, therefore suffice it to know the expression of in the special 
limit X — y — xq. We start from the analytic expression of the NLO correction to the 
hard-scattering kernel tabulated in Ref. [39]. After some straightforward manipulations, 
and paying particular care to the spurious singularity affiliated with the limit x y, we 

^ A comprehensive diagram-by-diagram analysis for the NLO correction to n EM form factor can be found 
in Refs. [29, 39]. 
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end up with a relatively compact expression: 

^ii' ^. f . I) 4 2 .n . . .„ I) . ^^(3 . 2 in .) 1,> I . iu.(.) 

4, 2 1 - 32x + 157x2 , 1,2 4 - 13x, 1 - 102a; tt^ 
+ 3'"^+ 36;^ '"^+6'"^+^6^'"^ + ^6^-36- '^^^ 

which is asymmetric under the interchange x ^ x. Note this expression is valid specifically 
for Nc = 3, and we have not attempted to keep the full track of the general Nc dependence 
in all the occurring color factors. Obviously, the ln{Q^/iJ,jp) term in (22) is the trace of the 
coUinear IR singularities encountered in the original NLO quark amplitude. 



2. NLO correction to the jet function: 4>^^^ 

According to (17), an important class of NLO correction is from convoluting the tree-level 
hard-scattering kernel with one NLO Be distribution amplitude and one LO Be distribution 
amplitude. A novel feature in our light-cone treatment is that, the jet function associated 
with the Be LCDA can be systematically improved in perturbation theory. It is the very 
feature that renders the strategy of refactorization practically useful. 

Fortunately, the NLO perturbative correction to the jet-function for the Be meson, accu- 
rate at the LO in v expansion, is not needed to be computed in this work, because recently 
it has already been calculated by Bell and Feldmann [22] Here we just quote their result: 



2 V MUxo-xy 



Xo XXo \ Xo Xo 



XX \ 1 , \ f ^ - Xo 

+ - d{x-xo){ 2xoXo In — 



+ z^] +-S'{x-xo){2xoXoln^ + xo-xo)} . (23) 



Here the and -prescriptions are understood in the sense of distributions. For a 

test function f{x) which has smooth behavior near x — xq, its convolutions with the 



In Ref. [22], the authors carried out the calculation by taking K meson as the prototype for a nonrelativistic 
system composed by s and u quarks, with m„ ^ nig. They performed a separate study for the Be meson, 

assuming that it belongs to the heavy-light meson family, alike to . Since the Be meson is regarded as 
a truly nonrelativistic bound state in the current work, we thus transplant their expression of cj)^^^ for the 
K meson, rather than theirs for Br. 
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and "++" -functions are given by 

fdx [g{x)Uf{x) = Cdxg{x) {f{x) - f{xo)) , (24a) 
Jo Jo 

^ dx [g{x)]++f{x) = 2 dxg{x) (^f{x) - f{xo) - f{xo){x - Xo)^ . (24b) 

As one can readily tell from (23), the jet function (()^^\x) has a nonvanishing support in 
the full range from to 1, which is in sharp contrast to the infinitely-narrow LO jet function. 
In particular, the NLO jet function has developed a long tail, as a consequence of reshuffling 
the momentum fraction between c and 6 to a highly asymmetric configuration through ex- 
changing energetic coUinear gluon. It is interesting to compare this perturbatively-generated 
broad profile of the quarkonium jet function to the phenomenologically determined quarko- 
nium LCDA, which has a much narrower width of 0{v). 

As is well known, the evolution of the leading-twist LCDA of a meson is governed by the 
renormalization group equation, which is commonly referred to as Efremov-Radyushkin- 
Brodsky- Lepage (ERBL) equation [40, 41]. Particularly, the jet function of the Be meson 
obeys the following evolution equation: 



^ 4>(x,^il) rdyVo(x,y)^(y,^il) + 0(a% (25) 



2 ^ /-i 



din lip 
where 



Vo{x,y) = ^ 



1 + ]e{x-y) + -(l + ^— ]e{y-x) 



(26) 



1-y \ x-yj y \ y-x^ 

is the corresponding evolution kernel. One can explicitly check that, upon substituting (23) 
into (12), Eq. (25) is indeed satisfied. 

With the knowledge of the Born-order hard-scattering amplitude (18) and the LO ex- 
pression of the Be LCDA, ERBL equation can be utihzed to identify and resum the leading 
coUinear logarithms to all order in ag [19]. In particular, it is not difficult to apply ERBL 
equation to infer the single-collinear- logarithm scaling observed in Fig. 2. 

When convoluting the jet function with the hard-scattering kernel in the light-cone frame- 
work, one often encounters the inverse moment of x {x), i.e., the integral of the jet function 
weighted by 1/x (l/x), as indicated in (20). With the exphcit expression (23), it is straight- 
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forward to deduce the first inverse moments of x and x at NLO accuracy: 







X 



4xo 



(3 + 21nxo)ln-5|- + 4Li2(xo) 



Ml 



2 In^ xq — 2(1 + 2>Xq) In xq — Qxq In a;o + 6 



27r2 



{x 



da; 



(a; 



(27a) 
(27b) 



so 



The respective first inverse moments at LO accuracy are trivial, i.e., 

The first inverse moment of x at NLO depends on factorization scale logarithmically. It 
is worth noting that, the ln(/i|n/M^^) term has the same coefficient as the ln((5^//i|') term 
in the NLO hard-scattering kernel in (22). This is guaranteed by the general principle of 
coUinear factorization framework. 



3. NLO correction to Be decay constant: f^j 

As outhned in (20), the last missing piece for a complete NLO analysis is the NLO 
perturbative correction to the decay constant. This information can be inferred through 
matching the QCD axial vector current onto its NRQCD counterpart to NLO in ag, which 
has also been available long ago [42] : 

fS = -lcF + lc4xo-xo)ln^, (28) 

Z 4 Xq 

which is symmetric under xq Xq. 



4- Final NLO prediction to pjjl 

Now it is time to piece all the relevant elements together. Substituting 'X^^\xo, xq) given 
in (22), (x^-^)^^^ in (27), and f^] given in (28), into equation (20), we finally obtain the 
complete NLO light-cone prediction to the B^. EM form factor: 



p(l)/n2^ 27rCFQ;.(/i|) J ec 
-^Lclv ) - 







^f|-21na;o + ln^V%(3 + 21na;o)ln 



+ 3U(x)+^ln'x I 1-^2x0 + 37x1 ^ - 85xo 1 - 102xo 177r^ 

+ 6 ui2[Xo) + - \.n Xq-I ——^ mxoH mxoH — -r— 

D 6bxQ 6bxo 6bxo 6b 

- (Cc 66, Xo ^ Xo)> . (29) 
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Note that the factorization scale /jLf has disappeared from this ultimate NLO prediction, as 
it should be. This happens because the IR logarithm \n{Q^ / nj^) in smoothly merges 
with the UV logarithm ln(/x|,/Mg ) in {x~^)^^\ The simple manifestation of the coUinear 
logarithms highlights one attractive power of the light-cone approach. 

Equation (29) constitutes the climax of this work. This NLO prediction from the refac- 
torization approach is impressively compact, in stark contrast with the extremely involved 
NLO expressions from the NRQCD approach, though they are supposed to be completely 
equivalent at sufficiently large Q^. 

D. Comparison of light-cone and NRQCD predictions at NLO in Ug 

1. Physical Be 

We arc now in a position to make a critical comparison between the NRQCD and light- 
cone predictions to the B^. EM form factor through NLO in a^. 

In the left panel of Fig. 2, we have depicted the ratio -PlcV-^lc ^ ^ function of Q, which 
is generated according to (29) and (19). The values of xq and xq are fixed by taking rric — 1.5 
GeV and rrib — 4.8 GeV, relevant for a physical state. The corresponding NLO NRQCD 
prediction, has also been juxtaposed in the same plot. As anticipated, these two predictions 
do converge together as Q becomes much greater than Mb^- hi Q = 10^ GeV, the light-cone 
asymptote is about 9% lower than the NRQCD result; but at Q = 10"^ GeV, the relative 
error between the light-cone and NRQCD predictions shrinks as tiny as 3 x 10~^. 

This highly nontrivial agreement has several important implications. First, the predic- 
tion from the light-cone approach armed with the refactorization strategy, (29), can be 
unequivocally regarded as the asymptotic expression of the NLO NRQCD prediction. This 
identification is expected to hold order by order in ag. 

This agreement also corroborates, in an indisputable manner, the correctness of the NLO 
NRQCD calculations presented in this work, the correctness of the NLO calculation in the 
light-cone side, i.e., the hard-scattering kernel [39], the jet function for the Be meson [22], 
as well as the NLO correction to the Be decay constant [42] . 

We have stressed that, due to the pathological complexity of the analytic NLO expression 
from the NRQCD side, it would be extremely time-consuming, if not impossible, to deduce 



22 



the asymptotic behavior directly from this expression itself. Remarkably, by proceeding 
along a rather different route, we are able to deduce, with much ease, the desired asymptotic 
NRQCD expression in closed form. 

Lying in the heart of the utility of refactorization, is its maximal exploitation of scale 
separations. Hard exclusive reaction involving quarkonium is complicated by the coexistence 
of several widely-separated energy scales, e.g., Q, m, mv, and so on. Refactorization strives 
to dissect such a multi-scale problem into several simpler steps, each of which only focuses 
on a single scale and can thus be tackled easily. Prom the example of Be EM form factor, 
we hope to have convinced the readers that this strategy is indeed much more efficient than 
the brute-force calculation within the NRQCD factorization. 



2. Fictitious Be 

To sharpen our understanding, it is also instructive to examine the EM form factor of 
a fictitious Be meson. We have chosen rric — mi, — 3.15 GeV in such a case. In the right 
panel of Fig. 2, the ratio F^^/F^q for the fictitious Be state is also shown for a wide range 
of Q. As can be clearly seen, the light-cone prediction starts to overlap with the NRQCD 
prediction since Q > 10^ GeV. At Q = 10^ GcV, the light-cone asymptote is about 3% 
lower than the NRQCD result; but at Q = 10^ GeV, the fractional error between these two 
predictions reduces to 6 x 10~^. 

For a fictitious Be meson, the asymptote of the NLO NRQCD predictions has already 
been known analytically in Sec. Ill, thus we can also make a comparison at the analytic 
level. Substituting xq — xq = ^ into (20), the light-cone prediction becomes 

- "^'h' ( I i i, + {x-'r - 4. (30) 



where the last entity is due to f^] for the equal mass case. The intended expressions for T^J^ 
and {x~^)^^^ can be simply deduced from (22) and (27) The ratio of the NLO prediction 



Note that the jet function for the fictitious Be state is equal to that for r/c. The NLO perturbative 
correction to the jet function for r]c was first evaluated in [21]. However, as pointed out in Ref. [22], the 



(f>^^\x) determined in [21] does not respect the due normalization condition. 
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to the LO one for the fictitious Be state reads: 



4,2^ 25, ^ 25 2%^ 

--ln^2H In 2 , (31) 

3 9 9 9' ^ ^ 

where the artificial factorization scale fip cancels. This expression can also be directly 
obtained by substituting xq — ^ in (29). Inserting /3o — ^ (for i^-f — 5), this hght-cone 
prediction, reassuringly, agrees exactly with the asymptotic NRQCD result tabulated in (8). 

Prom the light-cone perspective, now it should be clear why the coefficient of the coUinear 
logarithm in the asymptotic NLO NRQCD expressions for the EM form factor of the ficti- 
tious Be and the rjc—'y form factor, appears to be proportional to Cp (3 — 2 In 2). This occurs 
because the coUinear logarithms in both cases can be identified from the inverse moment 



E. Improving the scale dependence by refactorization strategy 

In the preceding analysis, our primary goal is to verify that, for the leading-twist hard 
exclusive process involving two quarkonia, exemplified by Be form factor, the light-cone 
approach, when armed with the machinery of refactorization, can be utilized as an efficient 
and elegant tool to systematically reproduce the asymptotic NRQCD prediction. As we 
have seen, our approach has withstood nontrivial test at the NLO level. We anticipate that 
the refactorization method works presumably to any fixed order in as- 

In higher-order calculation from the NRQCD factorization approach, it is a common 
practice to attach all the occurrences of with a single scale, fj,R by default. It has been 
observed that, for the NRQCD prediction of the double charmonium production process 
e+e" J/ip + r)c, the scale dependence does not improve at all even after including the 
NLO correction [11, 12]. This may be viewed as a serious drawback of the conventional 
NRQCD factorization approach, presumably attributed to the inadequate disentanglement 
of the scales Q and m in NRQCD short-distance coefficients. 

In the light-cone approach that implements refactorization, the complete NLO result 
encompasses contributions from several different ingredients. As stressed in Sec. IV A, the 
different as associated with the different source in (17), should in principle be evaluated 
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at different scales, e.g., /ip, H^r, and M^^, respectively. While it might be common to set 
A*ii = ~ Q, it is certainly more reasonable to evaluate the ag associated with the last 
piece in (17) at a lower scale around Mb^. By this way, one could, presumably, make a 
physically more sensible NLO prediction than the standard NLO NRQCD prediction. 

Since our light-cone approach has done a finer job in disentangling the distinct scales 
than the standard NRQCD factorization approach, it is natural to envisage that the scale- 
dependence might be significantly reduced with the implementation of the refactorization 
strategy. We hope to explicitly illustrate this attractive property in future publication. 

V. ASYMPTOTIC BEHAVIOR OF HEAVY-LIGHT MESON EM FORM FACTOR 

Thus far, the meson has been treated as a genuine quarkonium state. However, it 
has also occasionally been classified as a heavy-light meson by some authors. Consequently, 
we cannot resist the temptation to apply the NLO analysis performed in previous sections 
to explore certain features about the hard exclusive processes involving heavy-light meson. 
For definiteness, in below we will take the asymptotic behavior of the meson form factor 
as a concrete example. This section has somewhat digressed from the main thread of this 
work, hence an uninterested reader may skip this and go on to next section. 

Obviously, NRQCD factorization approach, which is tailor-made to describe the hard 
processes involving quarkonium, will no longer be a rigorous approach to deal with hard 
exclusive process involving B meson. Applying NRQCD factorization approach literally to 
analyze the B meson form factor amounts to modeling the B meson as a nonrelativistic 
bound state composed of a heavy h quark and a light constitute u quark 

Provided that is asymptotically large, the light-cone framework should still apply, 
but the strategy of refactorizing the B meson LCDA as the product of the B meson de- 
cay constant and a perturbatively calculable jet function looks obviously unjustified, since 
the necessary condition mb, ^ Aqcd has been violated. On the other hand, although 
modeling the B meson LCDA as ^b{x) — ^^§^ ^{^ ~ ^o) may seem to be a overly rough 



approximation, it may arguably capture some reasonable physics provided that m„ is taken 
as the constituent quark mass of order Aqcd- That is, this approximation is compatible 

For a phenomenological investigation of time-like D^*^~^D^*^~ form factor from this perspective, see [43]. 
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with the characteristic picture of a heavy-hght meson, that the momentum carried by the 
hght quark is predominantly soft. 

This said, let us recklessly apply the NRQCD factorization and our refactorization ap- 
proach to the B meson form factor at large Q^, bearing in mind that a constitute quark 
model picture has been assumed. First observe that in the heavy quark limit m„/mf, 
(hence — >■ and xq 1), the diagrams in which the EM current is attached to the up 
quark line can be neglected with respect to those in which the EM current coupled to the b 
quark, due to the much higher virtuality of the exchanged gluon in the former case. Starting 
from either (6) or (19), one can readily infer the asymptotic expression of the LO S+ form 
factor: 



where fs is the S+ meson decay constant. Besides the normal scahng, there is a 

pronounced enhancement factor brought in by I/xq. Obviously, an extra complication of this 
case with respect to the B^. form factor, is that there emerges one additional nonperturbative 
parameter, 1/xq. 

At first sight, one may view (32) as a naive estimate from the nonrelativistic constitute 
quark model, which should not be attached with too much significance. 

One alternative formalism, yet tailor-made to deal with hard exclusive reaction involving 
heavy-hght meson, has been developed by Grozin and Neubert some time ago [44]. They 
introduced a pair of new B meson nonperturbative distribution amplitudes 0±(c<;), specifi- 
cally defined in the context of heavy quark effective theory (HQET), where u denotes the 
light-cone energy of the spectator quark. A nonperturbative parameter that appears in vir- 
tually every exclusive B meson decay process, is the first inverse moment of the B meson 
"leading-power" distribution amplitude the so-called 1/Xb = / duj(j)+{uj)/uj, which scales 
as 1/Aqcd by dimensionless counting. 

The formalism developed in [44] (see also [45]) has later found interesting apphcations 
in the high-energy B and D meson production processes [46]. Analogous to the NRQCD 
factorization being a proper framework to describe quarkonium-involved hard reactions, we 

The term power correction in this context refers to effects suppressed by powers of 1/mb, which should 
be distinguished with those by powers of 1/Q. 




(32) 
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FIG. 3: The ratio F^^l (Q^) / F^^l (Q^) as a function of Q with Mb+ = 5.28 GeV, nf=5 (/3o = f ), 
and fiji = Q. We have chosen m„ = 0.48 GeV and mf, = 4.8 GeV. The NRQCD factorization 
prediction is represented by the sohd hne, while the hght-cone predictions from (29) by the dashed 
hne. The dot-dashed hne is depicted according to (33). Numerically, the dashed line can be 
parameterized by — 2. 394 In Q + 9.096, and the dot-dashed line by — 2. 394 In Q -|- 8.393, where Q is 
in the unit of GeV. 

may promote the theoretical framework underlying [44] as the HQET factorization approach, 
which may be suitable to tackle the heavy-light-meson-involved hard reactions. 

As the NRQCD factorization is based on v expansion, the appropriate expansion param- 
eter in HQET factorization is Aqcd/''^6- In [46], a useful connection between the NRQCD 
model calculation and the HQET factorization approach was offered: for a heavy-light me- 
son production process, first obtaining the reaction amplitude by employing the NRQCD 
factorization method, then substituting the singular factor I/xq by rrib/XB while neglecting 
all the remaining occurrences of Xq elsewhere. A gratifying fact is that, once such a recipe is 
adopted in (32), one then correctly reproduces the asymptotic LO prediction made in [44]. 

By far, a NLO analysis to the B meson form factor based on the HQET factorization 
framework has not emerged yet. Hence it may be interesting to examine the asymptotic 
behavior of the B meson form factor at NLO in our NRQCD factorization model and light- 
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cone model. 

In Fig. 3, we have shown the NRQCD model prediction and the corresponding hght- 
cone prediction to the ratio of F^^\Q^) / (Q^) over a wide range of Q. As expected, 
one observes good agreement between these two approaches after Q > 200 GcV. A distinct 
feature of this ratio is that it decreases as Q increases, which is opposite to the case of the 
Be form factor. 

To see such a scaling behavior in a more lucid way, one may further utilize the heavy 
quark limit, m„/mb — > 0, to simplify the light-cone prediction (29). After some algebra, we 



One can see from Fig. 3 that such a limiting behavior is numerically close to the exact 
NRQCD and hght-cone predictions at large Q'^. From this equation, it can be easily 
understood why the collinear logarithm now has developed a negative slope, because of 
mu/mb < 1. 

Although the prediction (33) is far from being rigorous, it may still contain some essen- 
tially relevant ingredient. For example, the coefficient of the collinear logarithm depends 
on Iuxq. Conceivably, this is what one would expect from a NLO calculation in the HQET 
factorization, once one identify this term with the logarithmic moment of the (f)^{u), the 
quantity called cxb [47]. 

From our lesson of bridging the NRQCD and light-cone approaches to effectively describe 
the quarkonium production, it may sound appealing to ask whether a similar refactorization 
strategy can also be apphed to the heavy-light meson production. That is, is it possible to tie 
the HQET factorization and the hght-cone approaches fruitfully? Since the B meson LCDA 
still contains the collinear degrees of freedom of higher virtuality of order m^, it is natural to 
expect that such short-distance effect can be separated from the remaining nonperturbative 
part. If this reasoning works, the leading-twist B meson LCDA may be, conceivably, factored 
into the convolution of the perturbatively calculable jet function with the (p+{uj). One 
attractive point is that, such a refactorization program enables one to manifestly disentangle 
the perturbative collinear logarithm of ln{Q/MB) from the nonperturbative soft logarithm 
of ln(MB/AQCD)- While the former logarithm is controlled by the ERBL equation, one can 
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invoke Lange-Neubert evolution equation [48] to deal with the latter. 



VI. CHALLENGE OF THE REFACTORIZATION PROGRAM TO J/il) - rjc EM 
FORM FACTOR 

In recent years, one of the most widely-studied hard exclusive reactions is perhaps the 
double charmonium production process e+e" J/ip + rjc, which was observed at the B 
factories a number of years ago [5, 6]. This process provides a powerful probe to extract the 
J/ijj + r]c EM form factor in the time- like region. One can define this form factor, which be 
referred to as G{Q^) hereafter, as follows: 

W(P, e(A)) + VciP'MM = i G(Q') e^^^'^ P,Plel{\) , (34) 

where Q — P + P', and e(A) represents the polarization vector for the J/ip with helicity 
A. This specific Lorentz structure is constrained by the Lorentz and parity invariance, from 
which one can easily see that the outgoing J/ip must be transversely polarized, i.e., A = ±1. 

In analogy with (2), according to the NRQCD factorization, one can express the form 
factor GiQ"^) as 

G™«co(Q^) ^ „j iim^pM . ^ . (35) 

\J2Ncmc \J2Ncmc 

where we have only retained the contribution at LO in v. 

The NRQCD short-distance coefficient, C{Q\ rric), can again be organized in power series 

of the strong coupling constant, C = C*^"^ + ^C*^"*^-' + • • ■ . There are four lowest-order 
Feynman diagrams for hard-scattering process, and the tree-level result of C^^^ is [8] 

C^'\Q;m,) = 2567reeC^a,(/i^) (36) 

As is well known, since this process necessarily violates the helicity selection rule [49], so 
the form factor has to be suppressed by extra factors of 1/Q relative to the Be form factor. 
From (36) one can infer that the helicity conservation is violated by the heavy quark mass. 

The NLO perturbative correction to this process was first computed in NRQCD factor- 
ization context in [11], in which only the numerical result a,t Q = ^/s — 10.58 GeV of 
S-factory energy was presented. This calculation was later redone by Gong and Wang, as 
the maiden application of the FDC package [12]. They confirmed the result given in [11], 
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lO' 10^ lO' 

Q (GeV) 

FIG. 4: The ratio of time-like J/V- + ??c EM form factor G^^'> (g2)/G(o) (Q^) as a function of Q. We 
have adopted rric = 1.5 GeV, nf=A (/3o = and fiR = Q. Both of the exact and asymptotic 
NLO results in NRQCD factorization framework are shown, where the former is represented by 
the solid line, and the latter [defined in (37)] by the dashed line. The vertical line marks the place 
at Q = y/s = 10.58 GeV of B factory energy. Numerically, the dashed line can be parameterized 
by 2.1671n2 Q - 3.077 In Q + 6.672, where Q is in the unit of GeV. 

in addition, they also provided the fully analytical expression for C^^\ This information is 
certainly useful, since one can then acquire a global understanding of this form factor at 
asymptotically large Q^. After some judicious but tedious manipulations on the analytic 
expressions assembled in [12], we can extract the asymptotic behavior of this form factor at 
NLO in as'. 

MC^UQ)] ^ l^M.^_ll(21n2-l)ln^ + ^ln4 
dSm(Q) 24 m2 24^ ' ml 4 

71 59 , 2 23 TT^ 

+ — ln2 + — ln^2 . (37) 

8 24 18 36 ^ ' 

where we have normalized C*-^^ with respect to C*-"-*, and suppressed the imaginary part 
for simplicity. In Fig. 4, we have depicted the ratio of the time-like J/ip + r]c EM form 
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factors G'(i)(g2)/G'(o)(Q2)^ or equivalently, C^^\Q'^) / C^^^Q"^) , as a function of Q, including 
both the exact and asymptotic NLO results. It is interesting to observe that, unlike in 
the Be case, the asymptotic expression starts to decently reproduce the exact NLO result 
at rather low scale of Q, that is, the mutual agreement at 5- factory energy seems already 
quite satisfactory 

A peculiar feature is that, as manifested in Fig. 4, here the leading scaling behavior of the 
NLO correction is actually governed by the double logarithm of type {Q^ / ml) This 
is in stark contrast to the asymptotic behavior of the Be form factor, as given in (29) and 
(31), whose leading behavior at NLO in is represented by the single coUinear logarithm 
only. 

It is interesting to investigate the numerical significance of this double logarithm term. 
Setting jiR — Q and rric =1.5 GeV in (37), we find the value of the double logarithm ^ In^ £2 
equals 8.27 aX Q — 10.58 GeV, which already constitutes 72% of the full asymptotic result 
ReCi^ym/Ci^ym = 11.47. The dominance of this double logarithm becomes very prominent 
at higher energy. For example, if the e'^e" center-of-mass energy is chosen to be near the 
pole, i.e., Q = 91.19 GeV, one finds that the double logarithm term becomes 36.55, 
which is 99% of the the full asymptotic result ReCisym/Cisym = 36.91. This numerical study 
undoubtedly indicates that, the importance the double logarithm term may severely ruin 
the stability of the fixed-order perturbative expansion, and in order to obtain the controlled 
prediction, it seems compulsory to resum these types of logarithms in NRQCD short-distance 
coefficient to all orders in ctg, which is perhaps relevant even at B factory energy. 

At first sight, the occurrence of this double logarithm seems to be at odds with the general 
principle of light-cone framework, especially our refactorization program. Recall that when 
computing the 0{as) quark amplitude for tt form factor, the double IR poles and the single 
IR soft poles have canceled away upon summing all the NLO diagrams, and only the single 
coUinear (mass) singularities can survive in the final answer. 

In passing, we note that the rising of double logarithm are not something completely 

To be concrete, at Q = 10.58 GeV, one finds that ReCisym/Casym = 11.47, which is only 3% larger than 
the exact NLO result 11.15. 

It is interesting to note that Ini[ci,sym(Q)]/ci,sym(Q), whose full analytic expression is not given here, 
contains the single logarithm —tt^Ih^. Conceivably, with this input, one can reconstruct this double 
logarithm in the real part in (37) via the dispersion relation. 



31 



new for the double quarkonium production processes in NRQCD factorization approach. 
In particular, this type of scaling behavior has already been encountered in the exclusive 
processes such as T — >• J/il) + 7]c [50] and 77^ — )■ J/il) + J/'^ [27, 51]. 

One can quickly see that, all the aforementioned doublc-quarkonium reactions share a 
common trait, i.e., that the celebrated helicity selection rule has been violated in all of 
them. This implies that the asymptotic scaling behaviors of these reactions are suppressed 
by powers of l/Q with respect to that of the Be form factor. Thus to make a leading 
nonvanishing prediction for such type of processes in light-cone formalism, one necessarily 
needs to include the higher-twist LCDAs of quarkonium. 

Unfortunately, a long-standing difficulty associated with these higher-twist LCDAs in 
collinear factorization framework is that, one will encounter the ubiquitous end-point singu- 
larity when convoluting these LCDAs with the hard kernel [2] . By far there is no universally 
accepted recipe to remedy this notorious problem other than some ad hoc phenomenological 
parametrization For example, there seems no consistent way to investigate the NLO 
correction for the fm EM form factor in the light-cone formalism. By the same token, this 
may cast some shadows on the solidity of those phenomenological studies on j* ^ J/ip-\- rjc 
from the light-cone approach [13]. 

We are thus facing a dilemma to carry out the refactorization procedure to the processes 
like 7* — >■ J/V' + rjc, because of our incapabihty of making a consistent NLO analysis of this 
higher-twist reaction in the light-cone framework 

Notwithstanding the theoretical hmitation of the light-cone approach, there is still some- 
thing worth learning solely based upon the NLO NRQCD results. After a diagram-by- 
diagram anatomy of the FDC output [12], one is able to single out those diagrams that 
contain double logarithms, which are depicted in Fig. 5. It is easy to observe that the sub- 

The so-called zero-bin subtraction method [52] has been recently proposed to solve the endpoint singularity 
problem. But there seems to exist some controversy about its validity. 
^"^ Recently there has come out an all-order-in-as factorization proof for the processes of e+e~ annihilation 
into double quarkonia [53], which seems to largely build upon the arsenal of collinear factorization. The 
key technique in their proof is the soft and collinear approximations of QCD interaction which conserve 
the quark helicity. Hence their proof may be regarded as applicable to those leading-twist (helicity- 
conserving) reactions such as Be EM form factor. For higher-twist (helicity-suppressed) processes such as 
e+e~ J/tjj -\- r]c, one has to include the spin-dependent interactions to flip the quark helicity, which is 
necessarily beyond the soft and collinear approximations employed in [53] . Therefore it is not clear to us 
whether their proof can be safely applied to e+e~ — >■ J/ip -\- r/c or not. 
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f) 
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J) 



OOOOOO' 



5) 
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FIG. 5: The NLO diagrams that contain the double logarithm ln^((5^/m^) (in Feynman gauge) 
for the process 7* — )• J/ip + rjc- The respective charge-conjugated diagrams have been omitted. 



set of ladder diagrams [Fig. 5a) through c)] has a color structure oc ii{T'^T'^T^T^) oc C|n, 
while the remaining subset of diagrams has a different color structure oc tr(T'*T''T"T^) oc 
{Cp — \Ca)Cf, where Ca = Nc is the Casmir for the adjoint representation of SU (Nc). One 
can examine that the sum of all the diagrams in each subset still contains a nonvanishing 
double logarithm. As a result, for an arbitrary Nc, the coefficient of the double logarithm ap- 
pearing in (37) is a linear combination between Cp and CfCa- This color-factor dependence 
differs from the standard Sudakov double logarithm. 

It is instructive to reexamine the same topologies of diagrams in Fig. 5 for the case of 
vr form factor. As is illustrated in [29, 39], in that case the diagrams Fig. 5a) through c) 
only contain a single IR pole, in contrast to the case of 7* — )■ J/ip + rjc- Furthermore, for 
the remaining diagrams in Fig. 5, even though each of them still contains double IR pole, 
their sum does not. This comparative study clearly shows that the double logarithms are 
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sensitive to the helicity structure of the specific process. 

It is perhaps more transparent to trace the origin of these double logarithms by employing 
the method of region [54]. Following the version of refactorization outlined in Ref. [19], one 
may identify the relevant degrees of freedom in the NRQCD short-distance coefficient C^^\ 
There is always a ^^hard" region, in which the loop momentum scales as p'^ ~ Q. In this 
region, the mass of heavy quark should be treated as a small perturbation. There are several 
kinds of "infrared" modes, soft {p^ ~ m), collinear (p"^ ^ Q,p~ ^ rri^/QiPi. ~ rn), and 
anti-collinear (p"^ ~ /Q,p~ ^ Q,p± m). The mass of the c quark must be retained 
in these lower-energy regions. The validity of NRQCD factorization guarantees there is no 
overlap between these "infrared" quanta and those truly infrared modes intrinsic to NRQCD 
such as the potential mode. In the case of Be form factor, the leading contributions only 
arise from the hard and collinear (or anti-collinear, but not simultaneously both) regions, 
while the soft region cannot make a net contribution at leading power. In the J/ip + rjc 
case, it is expected that each of these four regions will make a contribution in the leading 
nonvanishing power (of course, there should be an overall power-suppressed factor). The 
double logarithm should originate from the overlap between the (anti-)collinear and soft 
regions. In this sense, we may dub this type of double logarithm as the power-suppressed 
Sudakov logarithm. 

As has been expounded in [55], wherever the double logarithm appears, the respective 
loop integrals in the two overlapping regions become ill-behaved separately. Rather one needs 
to introduce the extra regulator besides the dimensional regularization for each individual 
region, such as the analytic regularization. This artificial regulator will be eliminated only 
after summing the contributions of both regions. This symptom signals the breakdown of 
the naive coUinear-soft factorization, which is the very cause for the endpoint singularity. 
In phenomenological study of the p7r(or J/iprjc) EM form factor in light-cone approach at 
LO, one may dehberately choose the profile of the LCDAs such that it falls off sufficiently 
fast near the endpoint, as an means to circumvent the end-point singularity problem (for 
example, see [13, 56]). However, the real problem is that, this symptom is deeply rooted 
in the breakdown of the collinear factorization for higher-twist reaction, which cannot be 
simply overcome by invoking some phenomenological trick. 

Even though it is feasible to reproduce the asymptotic behavior in (37) with the aid of the 
method of region, one perhaps still can not proceed far without a general guidance of a valid 
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factorization theorem. For instance, it still remains to be a harsh challenge to identify and 
resum these power-suppressed Sudakov logarithms to all orders in as- Perhaps a satisfactory 
answer to this question demands that one first resolves the end-point singularity problem in 
a consistent manner. This direction definitely deserves further exploration. 

VII. SUMMARY 

It is an indisputable fact that both the NRQCD factorization and the coUinear factor- 
ization approaches have their own strengths and limitations in describing hard exclusive 
reactions involving heavy quarkonium. In this work, we have illustrated how to coherently 
tie these two approaches to achieve the optimized predictive power. This is made possible 
by invoking the strategy of refactorizaiion, e.g., by further factoring the quarkonium LCDA 
into a sum of the products of the universal yet pcrturbatively-calculable jet functions and the 
nonperturbative vaucuum-to-quarkonium NRQCD matrix elements (i.e., quarkonium decay 
constant). Through a comprehensive comparative study, we have verified that, for a class of 
hard exclusive process involving two quarkonium, exemplified by the Be form factor at large 
Q^, the light-cone approach with refactorization can be utilized to reproduce the (leading- 
twist) asymptotic result of NRQCD factorization prediction, at the NLO in while at the 
LO in V. We hope to have convinced the readers that this refactorization program is much 
simpler, and, more efficient, than the brute-force NRQCD-factorization-based calculation. 

Quite conceivably, for most realistic quarkonium production processes, e.g., rjc — j tran- 
sition form factor to be the simplest, it is beyond our current technical capability to ana- 
lytically investigate the next-to-next-to-leading (NNLO) perturbative correction in NRQCD 
factorization formalism. Fortunately, the refactorization approach, owing to its maximally 
disentanglement of the scales, will serve as the indispensable, and, perhaps the only viable, 
calculational device to fulfill this goal. Beside such a remarkable technical advantage, the 
strategy of refactorization also helps one to gain better theoretical control, e.g., to resum 
coUinear logarithms to all orders in as, and to significantly reduce the scale dependence 
for a fixed order calculation. We hope to explicitly show these attractive features in future 
publication. 

There are quite a few leading-twist quarkonium production processes for which the refac- 
torization program may be of interesting applications, e.g., B meson exclusive decay to 
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a S-wave charmonium plus a light meson, and 77 — > + B~ . It is also worth further 
factorizing the quarkonium LCDA to higher order in v expansion, by which one is then ca- 
pable of elegantly reproducing the NRQCD predictions for P-wave quarkonium production 
or relativistic correction to S'-wave quarkonium production. 

Another interesting direction is to pursue a similar refactorization program in hard ex- 
clusive reactions involving a heavy-light meson, e.g., to further factorize the leading- twist B 
meson LCDA into the convolution of a perturbatively-calculable yet universal jet function 
with the B meson distribution amplitude specifically defined within the HQET context. 

It is an empirical fact that our refactorization program confronts serious obstacle when 
applied to helicity-suppressed double-quarkonium production processes, such as 7* J/'ip + 
Vc Vb ^ J li' + J li'i and T — 7> J/ip + 1]^. Wc believe that the double logarithm appearing in 
the NLO NRQCD short-distance coefficients in these processes is intimately linked with the 
end-point singularity problem in light-cone approach. In our opinion, a satisfactory control 
of this double-logarithm for these processes will offer a pivotal insight into the ultimate 
solution to this long-standing problem in coUinear factorization. 
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